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Nonlinear oscillations of a bubble carrying a constant charge and suspended in a fluid, undergoing 
periodic forcing due to incident ultrasound are studied. The system exhibits period-doubling route 
to chaos and the presence of charge has the effect of advancing these bifurcations. The minimum 
magnitude of the charge Qmin above which the bubble's radial oscillations can occur above a 
certain velocity ci is found to be related by a simple power law to the driving frequency uj of the 
acoustic wave. We find the existence of a critical frequency ujh above which uncharged bubbles 
necessarily have to oscillate at velocities below ci. We further find that this critical frequency 
crucially depends upon the amplitude Pg of the driving acoustic pressure wave. The temperature 
of the gas within the bubble is calculated. A critical value Ptr of Ps equalling the upper transient 
threshold pressure demarcates two distinct regions of uj dependence of the maximal radial bubble 
velocity Vmax and maximal internal temperature Tmax- Above this pressure, Tmax and Vmax 
decrease with increasing cj while below Ptr^ they increase with uo. The dynamical effects of the 
charge and of the driving pressure and frequency of ultrasound on the bubble are discussed. 

PACS numbers: 05.45.-a, 05.90.+m, 43.25.Yw, 43. 35. HI 



I. INTRODUCTION 

The stability and oscillations of a gas bubble suspended 
in a liquid under the influence of an acoustic driving pres- 
sure field in the ultrasonic frequency range have been 
the subject of a large volume of scientific literature [l|- 
ll4J . Studies on the system have been made from different 
viewpoints coming from its diverse applications and oc- 
currences. Ultrasound is routinely used in medical ultra- 
sonography including echocardiography, lithotripsy, pha- 
coemulsification, use in treatment of cancer and for den- 
tal cleansing. Other significant applications of ultrasonic 
forcing of fluids in which studies of bubble dynamics and 
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cavitation become very important are in sonochemistry, 
sonoluminescence, ultrasonic cleaning of materials, waste 
water treatment and in focussed energy weapons. Cavita- 
tion events which involve violent collapse of micron-sized 
bubbles in the fluid can cause immense damage to the 
surfaces they are in contact with. Studies of cavitation 
events in pumps, turbines, surfaces exposed to hydrody- 
namic flow, etc., continue to be of immense interest in 
industries and in technological designs of devices. 

Rayleigh's study of bubble cavitation was motivated 
by the need to understand and explain the damage 
to ships' propellers [1]. Under ultrasonic forcing, the 
behaviour of a bubble in a fluid depends heavily upon 
its ambient radius and the amplitude and frequency 
of the driving sound field. Thus the bubble can show 
regular oscillatory behaviour which can be periodic or it 



can show highly irregular oscillations which are chaotic 
and of unpredictable amplitude. For applications where 
damage caused on surfaces due to bubble cavitation 
can be disastrous, such as in medicine, it is desirous 
to operate the sonic device in a "safe" regime, and / 
or to be able to have control over the bubble's motion. 
Often in biological systems, it is known that bubbles in 
fluids can be electrostatically charged. Studies of the 
dynamics governing the oscillations, growth and collapse 
of charged bubbles are therefore of immense relevance 
because of their prevalence in diverse applications and 
situations. Experimental and theoretical work on the 
presence of charge on gas bubbles in fluids goes back 
to, for example, the work of McTaggart, Alty and 
Akulichev [ol, llO|, ll6|, ll7| , and more recently the work of 



ul, 



12|,Q. 



None of 



Shiran and Watmough and Atchley 

the work, though, has addressed the issue of dynamics 

of a charged bubble under ultrasonic forcing. 

It is interesting to know what effect the presence of 
electric charge on the bubble would have and see if the 
motion of such a charged bubble forced by ultrasound 
would vary significantly from that of an electrically 
neutral bubble in a fluid. This especially becomes of 
practical significance when we are looking at cavitation 
phenomena in fluids in real-life, be it in the context of 
cavitation in mechanical systems or in the case of bub- 
bles in fluids in living tissue in a medical context. Apart 
from the work in [15|, we are not aware of any other 
studies in the literature of the dynamics of acoustically 
forced charged bubbles suspended in a fluid. Their work 
however used the value 4/3 for the polytropic constant 
which entailed cancellation of all the charged terms; thus 
their work does not really address the issue of charge 
which it sets out to do. The extremely nonlinear nature 
of the system, and the presence of a large number of 
parameters do not facilitate a straightforward analysis 
and it becomes essential to take the aid of numerical 
methods to get an understanding of the dynamics 
governing the observed behaviour. In this work we 



report some studies on the dynamics of a charged bubble 
in a liquid (which we take to be water) when ultrasound 
is incident on it. We assume that heat transfer across 
the bubble takes place adiabatically, and the gas is a 
monatomic ideal gas. We therefore take the polytropic 
constant T = 5/3. 

In Section II we discuss briefly the nature of the radial 
dynamics of a charged bubble. Starting with a modified 
Rayleigh-Plesset equation, we obtain the time series of 
the bubble radius as also of its radial velocity and tem- 
perature. We also calculate the phase portrait of the 
bubble, under different pressure regimes. 
In Section III we discuss the pressure thresholds that 
influence bubble dynamics; we introduce the expansion- 
contraction ratio C which we had introduced in 23| that 
enables us to locate the presence of the Blake and up- 
per transient threshold pressures easily when plotted as 
a function of the driving pressure amplitude Pg . The ef- 
fects of driving frequency uj and charge on ( are demon- 
strated in the present work. 

The influence of Pg and uj on the bubble dynamics are 
investigated in detail in Section IV. We obtain an expres- 
sion for the minimum charge required on a bubble for ra- 
dial oscillations to occur at some velocity ci, as also the 
dependence on the forcing pressure amplitude Pg of the 
maximum forcing frequency ujh at which an uncharged 
bubble will oscillate with velocity ci. 
We then obtain, in Section V, the bifurcation diagrams 
for the system with driving frequency as the control pa- 
rameter, and also the bifurcation diagram with charge as 
the control parameter. We observe that the presence of 
charge on the bubble advances period-doubhng bifurca- 
tions with driving frequency as control parameter. In- 
creasing Pg causes the advancement of period doubling 
and halving bifurcations for charged as well as uncharged 
bubbles, and bands of chaotic behavior are observed at 
large Pg. 
The effect of charge and driving frequency on the max- 



imal temperature are discussed in Section VI. We note 
that the pressure regime in which the bubble is being 
forced (whether Pg is above or below the upper transient 
threshold pressure) determines the frequency dependence 
of the temperature, and we obtain rough limits on the 
maximum charge a bubble may carry depending on its 
ambient radius. 

We conclude the paper with a summary of the results in 
Section VII. 



the vapour pressure of the gas. We denote by P^ and 
uj = 27riy respectively {u being the driving frequency), the 
amplitude and angular frequency of the ultrasound forc- 
ing field. We consider water to be the liquid surrounding 
the bubble, and having density p = 998/c^/m^, viscosity 
77 = 10~^ Ns/m^^ surface tension a = 0.0725N/m^ 
and the velocity of sound in the liquid c = 1500m/s, 
Pq = 101/cPa, r is the polytropic index and e = 85eo, 
where eo is the vacuum permittivity. 



II. RADIAL DYNAMICS OF THE CHARGED 
BUBBLE 

In real-life situations, bubbles in fiuids often have some 
electric charge sticking to them. This has been seen in 
the case of gas bubbles in various liquids as well as for 
cavitation events in water. In our work, we adapt the 
procedures for describing cavitation and forced bubble 
oscillations (that has a long and extensive literature), to 
include the presence of charge. 

Description of ultrasonically forced bubble motion in a 
fiuid has been made through the Rayleigh-Plesset equa- 



tion 



m 



and its variants 



Ka yleign- rles 



modified 



to take into account compressibility of the fiuid or var- 
ious other factors. Proceeding as we did in our earlier 
work [23|, we further modify the form of the Rayleigh- 
Plesset equation for the evolution in time of the bubble 
radius R{t) employed by [19] to include the presence of 
a constant charge Q on the bubble as follows jl5| : 
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Rq denotes the ambient equilibrium radius of the 
bubble, Pq the static pressure, and Py = 2.34kPa, 



The modified Rayleigh-Plesset equation above can be 
simplified and rewritten in dimensionless form ^ 
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where r = R/Rq, t = ujt, P^y = P^/Pq, P^^s = Ps/Po 
and the overdot here corresponds to differentiation with 
respect to r, and where the following dimensionless con- 
stants have been used: 



c* = ^; F = 4^; H 
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We have employed the dimensionless form of the 
equations for obtaining their numerical solutions. In all 
the expressions that follow, and in the numerical results 
shown in graphs, we have rescaled the quantities by the 
appropriate factors and only displayed the dimensional 
form for a physical grasp of the magnitudes of the 
quantities involved. 

The presence of charge Q counters the effect of surface 
tension, reducing its effective value, and induces several 



interesting changes to the dynamics of bubble osciha- 
tions. In a previous work [23| we had obtained for the 
charged bubble, the Blake threshold and radius and 
also some results for the upper transient threshold for 
cavitation. In the following sections, we describe some 
interesting consequences of the presence of charge on a 
bubble. 

As the bubble expands and contracts, the surface 
charge density decreases or increases respectively. The 
presence of charge lowers the surface tension and for 
sub- micron sized bubbles, dominates over it, influencing 
the minimum and maximum values of the radius and 
the maximum velocities achieved by the bubble, and 
changing its point of collapse. A charged bubble achieves 
higher temperatures within it than an uncharged one, 
the collapse of the bubble being more violent in the 
charged case. 

The above results indicate that since the bubble os- 
cillations are more energetic for the charged bubble, the 
temperature attained by the gas within the bubble during 
its oscillations, would be higher as well. To confirm this, 
we calculate the temperature using equation (|3]) [2Q|. 

m = Tm{§^f\ (3) 

where h is the van der Waals hard core radius for the gas, 
h = R0/S.S6 for Argon [22]. This equation is obtained 
under the assumption that there is no exchange of heat 
from the gas to its surroundings, that the system is es- 
sentially adiabatic. 

This assumption is not strictly true as in reality the 
equation of state of the gas enclosed within the bubble 
can be either adiabatic or isothermal, depending upon 
the rate of collapse of the bubble and whether or not the 
various relaxation time-scales permit thermal diffusion to 
occur to and from the bubble. We use the expression in 
order to get an idea of the magnitudes achievable by the 
temperature in the presence of charge. To better visual- 
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FIG. 1: A. Time series of R, and B. Phase portrait (R vs R), 
for Ro = 5/im and bearing charge Q — 0(blue online) and Q — 
0.38pC (red online), Ps = (a) l.OPo, (b) 1.12Po, (c) 1.25Po. 
Amplitudes of R, R are larger; the phase portrait spans a 
larger space for the charged bubble than for the uncharged 
case (color online). 
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FIG. 2: Time series of temperature (in K) plotted as a func- 
tion of time for Rq = 5 /am at driving frequency u = 20kHz, 
and bearing charge Q = 0(blue online) and Q = 0.38pC (red 
online) for three values of Ps : (a) l.OPo, (b) 1.12Po, (c) 
1.25Po- Temperatures T are larger for the charged bubble 
than for the uncharged case (color online). 

ize the effect that the amplitude of the forcing pressure 
Ps and charge Q have on the bubble dynamics, we con- 
sider a bubble being driven at 20 kHz, i.e., the lower 
limit of the ultrasonic spectrum. Even in this lowest 
ultrasonic regime, the time series of bubble radius, ra- 
dial velocity, and temperature all show an enhancement 
in values due to charge. Moreover, Pg crucially deter- 
mines the dynamics of the bubble as illustrated in Fig- 
ures (jl|2I A-B. a,b,c). We have considered three values of 
Ps, Ps = l.OPo, 1.12Po and 1.25Po. These pressures are, 
respectively, below the Blake threshold Psiake^ at the 
upper transient pressure threshold Ptr, and above Ptr- 
As can be seen, the pressure regime in which the bubble 
dynamics occurs, crucially determines the behavior. At 
Ps = l.OPo, Tjnax ^ 370i^, the uncharged bubble tem- 
perature being marginally less than that for the charged 



bubble (Q = 0.38 pC); for P, = Ptr = 1.12Po, Tmax goes 
up to about 1520 K for the charged bubble and about 
1320 K for the uncharged case; and for Ps = 1.25Po, 
Tmax shoots up Still further, to about 24,000 K for the 
charged (and approximately 21,000 K for the uncharged) 
bubble. These temperatures vary by orders of magnitude 
and spell out the importance of Pg and Q. 

III. PRESSURE THRESHOLDS 

The Blake threshold determines the pressure threshold 
beyond which an acoustically forced bubble undergoes 
drastic expansion. After the Blake threshold and 
preceding the onset of bubble collapse following a larger 
threshold known as the upper transient threshold, Ptr, 
the bubble is essentially in an unstable regime. 
Depending upon whether the amplitude of the applied 
acoustic forcing pressure is greater or lesser than Ptr, the 
response of the bubble to the frequency of the applied 
pressure wave varies drastically. 

At low amplitudes of the forcing pressure (i.e., Ps < 
Ptr)^ increasing driving frequency causes a proportional 
increase in the bubble's maximum radial velocity Vmax- 
This happens upto some critical value of the frequency 
for that Ps after which Vmax rises more steeply but ac- 
companied with large oscillations. 

At larger amplitudes of the forcing pressure, with Pg ap- 
proaching the value of the Blake and upper transient 
threshold pressures, the situation is different. Vmax first 
decreases with increasing driving frequency upto a fre- 
quency uohci after which Vmax rises with frequency but 
with large oscillations. As could be expected from the 
above observations, a similar observation can be made 
regarding the maximum temperature Tmax of the gas in- 
side the bubble. 

A useful graphical illustration of the transient threshold 
pressures, i.e. of the Blake threshold {Psiake) and the up- 
per transient threshold {Ptr) pressures, can be obtained 



by plotting ( = {Rmax - Ro)/{Ro - Rmin) as a func- 
tion of the amplitude of driving pressure. This quantity 




FIG. 3: Plots of Rmin (left) and Rmax /Rmin (right) as func- 
tions of Ps for different values of ambient radius (for Q = 
and Q — 0.4 pC) at 20kHz forcing frequency, (color online). 

C, which we call the expansion-contraction ratio, handily 
shows the location of both the Blake and the upper tran- 
sient thresholds. Both these thresholds cannot be identi- 
fied easily at the same time from, for example, a plot of 
Rmax /Rmin as a fuuctiou of applied pressure amplitude. 
In Figure ([3|), the points of inflection of the curves corre- 
spond to the Blake threshold pressures for the respective 
Rq values. The effect of charge is clearly seen in reducing 
the threshold pressure as compared to the charged case. 
The upper transient threshold cannot be easily pinned 
down from this plot. While the Blake threshold is in- 
dicative of the threshold of the expansive growth of the 
bubble, the upper transient threshold demarcates where 
the violent collapse of the bubble occurs. A plot of C, ver- 
sus Ps for different values of Rq as shown in Figure (|4]) 
shows a rise of the curve till it peaks (at Ps = Psiake) fol- 
lowed by a trough or well (at Ps = Ptr) before rising up 
steeply for higher Ps (this has been discussed in some de- 
tail in our earlier work [23|). At pressures between Psiake 
and Ptr the bubble is in an unstable regime. This also 
explains the presence of large fluctuations or oscillations 
in the velocity vs. frequency plots at such intermediate 
pressures. The presence of charge shifts the threshold 
pressures to lower values. With increasing ambient ra- 
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FIG. 4: Plot of C = {Rmax - Ro)/{Ro - Rmzn) as a func- 
tion of forcing pressure amplitude Ps and Rq at 20kHz. The 
schematic at right shows the locations of the Blake threshold 
pressure and the upper transient threshold pressure on the 
C — Ps curve. Increasing bubble charge shifts the curve down 
and to the left, while increasing driving frequency decreases 
the steepness of the curve and flattens it. (color online). 

dius i?o, C loses its distinctive peak- valley appearance 
gradually. 

The maximum radius attainable by the bubble gradually 
increases with charge for a given driving frequency [23|. 
This can be understood from the fact that the presence 
of charge on the bubble decreases the effective surface 
tension. This causes the bubble to expand more easily 
in the negative pressure field. A casual reading might 
give rise to the observation that by the same argument, 
the minimum radius reached by the bubble would like- 
wise follow a similar trend, with Rmin for a charged bub- 
ble having a larger value than that of a neutral bubble. 
However, this is not so. It should be borne in mind that 
Rmin is influenced by the maximal velocity the bubble 
is able to reach. The greater the velocity, the smaller 
the Rmin that it collapses to. Hence, perhaps counter- 
intuitively, charged bubbles undergoing forced oscilla- 
tions, will achieve smaller values of Rmin than electrically 
neutral bubbles. 

Thus presence of charge leading to greater bubble expan- 
sion, in turn results in the bubble collapse being much 



more rapid and violent, shrinking the bubble volume 
more than in the case of the uncharged bubble. This can 
be seen in Figure (J3]) (left), where the minimum radius, 
Rmim reached by the bubble at the moment of collapse 
is plotted as a function of the driving pressure Pg and Rq 
for the charged and uncharged bubble. As was shown in 
greater detail in our earlier work J23| , Rmin reduces with 
increasing Q. 



IV. INFLUENCE OF AMPLITUDE AND 
FREQUENCY OF DRIVING PRESSURE FIELD 

The maximum radial velocity of the bubble attained 
during its collapse or contracting phase depends also on 
the driving frequency, the charge present on the bubble, 
as well as the amplitude of the driving pressure wave, as 
also on the initial radius Rq of the bubble in its quiescent 
state. Figure [5] are plots of the maximal radial velocity 
as functions of the driving frequency (a) and pressure Pg 
(b). There are several interesting features evident from 
the figures. 

The behaviour of Vmax above P^ = Ptr is different from 
that below it. The plots shown are for a bubble of Rq = 
5/im for which Ptr = 1.12Po. Fig.(j5^) shows that at Pg < 
Ptr^ Vmax iucrcascs as a function of driving frequency u 
while for Pg > Ptr Fig.dSb), it decreases. Increasing the 
driving frequency induces instability by producing large 
amplitude oscillations. 

For a given magnitude of pressure amplitude Pg, the 
magnitude of charge present influences the dynamical 
regime of the bubble. If the driving angular frequency of 
the applied pressure wave is a; at a certain pressure ampli- 
tude, for bubble oscillations to occur with some maximal 
radial velocity R = ci^ the charge present on the bubble 
should have some minimum magnitude Q = Qmini^)- 
At low frequencies, even an uncharged bubble might os- 
cillate at that velocity; however at higher frequencies, if 
charge Q < Qmin-, the radial bubble velocity would be 
smaller than ci. This is because as frequency increases. 
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FIG. 5: Maximum velocity vs (a) driving frequency, plotted 
at 7 different Ps values and for zero & non-zero values of 
charge; the Pg — 1.12Po curves (number 6, darkest [Q — 
black online) & just above that (Q — IpC, dark green online)) 
correspond to the upper transient threshold pressure Ptr (for 
Ro — 5/jm) and clearly demarcate distinct behaviour for Ps > 
Ptr and Ps < Ptr- (b) Vmax VS. Ps for different frequencies 
for Ro — 5/im, at 5 different driving frequencies (20 to 40 
kHz, increasing from left to right), for zero (•) and non-zero 
(Q = 0.4pC ,triangle) charge, (color online). 



the bubble does not get sufficient time to complete its 
expansion, so that its subsequent collapse occurs with 
smaller radial velocity than if expansion to a greater size 
had been done. The presence of charge reduces the sur- 
face tension and encourages expansion to larger radial 



dimension and the consecutive, more violent collapse to 
a smaller radius. 

That the change in Qmin{^) with uj show a bifurcation 
in the parameter space is clear from Figure [6l This tran- 
sition from zero to non-zero Qmin occurs at an angu- 
lar frequency ujh- For Pg = 1.35Po and Rq = 5/im, 
ujh = 23kHz. The magnitude of Qmin varies with driv- 

isovelocity curve 

(R: 




kinetic energy of the bubble would scale as the electro- 
static contribution ~ Q'^ / Rmin-, so that Q^ ~ Ci2ClRmin^ 

a2 being a prefactor with appropriate dimensions. Since 
for high applied pressures (that is, for values of Pg above 
the Blake threshold pressure or of the order of or above 
the upper transient threshold pressure) we know that the 
minimum bubble radius scales as the two-fifths power of 
the driving frequency Rmin oc o;^/^, it would follow that 
Q^ (X cio;^/^ so that Q ex cj^-^, which is close but not 
equal to the observed exponent of 0.25. Hence, this ar- 
gument, while it serves to give a lower bound for Qmim 
is insufficient. 

Proceeding more systematically therefore, we start by 
making a linearization of the Rayleigh-Plesset equation. 
Proceeding along the lines of [5] , the driving sound pres- 
sure is introduced through a small perturbation a, so 
that the total external field Pext can be written as: 



Po(l — acosuot) 



(5) 



The bubble oscillations R{t) about the equilibrium radius 
Rq can then be expressed as 



R = RQ{l^x{t)) 



(6) 



FIG. 6: Qmin versus frequency uj iso- velocity plots {v ~ 1500 
m/s) for a) Rq = 5iJ.m, for Ps = 1.35Po {Qmin is in pC, oj in 
kHz); b) schematic illustrating the nature of the iso- velocity 



plot, Qn 



[UJ-UJH 



^ ; c) plot showing the dependence of 



OJH on Ps for different values of Ro (curves 2, 3, 4, 5, 6, & 7 corre- 
spond to corresponding Rq values in /xm); d) plot illustrating 
variation of (dun/dPs) with Rq. (color online). 



ing frequency as 



Vr) 



a{uj -ujnf, 



(4) 



where the prefactor a has appropriate dimensions and 
depends on the value of the initial ambient bubble radius 
Rq, and h ^ 0.25. 

We could attempt to give a simple explanation for the 
frequency dependence of Qmin- We could argue that for 
a given value of constant maximal radial velocity ci , the 



where x{t) is a small quantity of order a. Substituting 
this equation in the Rayleigh-Plesset equation (1) and 
linearizing it, we get 



X -\- Px -\- UJqX 
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(7) 



where the damping coefficient /3, natural frequency of 
oscillation ujq of the bubble and F^xt are given by 
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Here only terms linear in x and its derivatives have been 
retained and (j)/ Rq is the equilibrium gas pressure in the 
bubble defined by 
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The particular situation of looking for conditions where 
the radial velocity is constant is thus implicitly satisfied. 
In eqn.©, Fext is given by 



F ^ — —P 

-^ ext — ^6 






X COS ujt — arctan 






[rI + if^) 



COS I Ujt — — 



2 c 



(10) 



where in arriving at the last line of eqn.([TO]), use has 
been made of the fact that ^^^ <C 1. Scaling the time 
as t = ujot for convenience, eqn.(7) can be solved exactly. 
Dropping the hat Q over t for convenience of notation in 
all of the following, the steady state part of the solution 
is found to be 
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where denotes the phase. 

Combining eqns.(ll) and (12) we obtain 
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(13) 
Again using eqn.(6) to rewrite x = R/Rq — 1 ,and 
V = R/Ro in eqn.(13), we obtain after some algebra an 
equation for uj: 
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After a careful look at each of the terms in this equation, 
we find that the dominant contribution of Q to duo /dQ 
occurs as a cubic term : 



duj 3 

dQ 



(17) 



as being a prefactor with appropriate dimensions. In- 
tegrating both sides of this equation between the limits 
corresponding to Q = and Q gives 



UJ - (jJh ^ asQ^ 



(18) 



where ujh is the frequency for the bubble with zero charge 
at which the Rmax = ci, so that 



/ / n1/4 

a^[uj -ujh) , 



(19) 



(the prefactor a'^ having appropriate dimensions), repro- 
ducing eqn.dlj that was obtained from an analysis of the 
numerical results shown in the plots in Figure (|6]). Hence 
it is very easy to predict the minimum charge Qmin re- 
quired on a bubble at a given applied pressure amplitude 
for attaining some particular value of the bubble's radial 
velocity, once uoh is known. 

There is another interesting feature to be noted in this 
transition. The value of the frequency uoh depends on 
the magnitude of the amplitude Pg of the driving pres- 
sure wave. Indeed, for a given ambient bubble radius i?o, 
uoh takes the simple linear form 



^H = hi + 62PS, 



(20) 
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FIG. 7: Q^ as function of Ps for a bubble of Ro — 2iim 
for three different values of driving frequencies v. The value 
of Ps at which Q/^ = 0, increases with increasing frequency. 
The physically reachable region is the area below the curve 
corresponding to Rmin > h. The region above the curve cor- 
responds to the (unphysical) regime where Rmin < h. (color 
online) . 

where hi and 62 vary with Rq. This can be seen clearly 
from the plot (Figure (J6])c). A further functional depen- 
dence of 62 on i?o, that is, of the slope dcon/dPs on Rq^ 
is also found (Figure (J6])d), and is of the form 



dP.. 



R^ 



-0.9 



(21) 



The maximal charge, Qmax-, which a bubble can carry, 
is bounded by the fact that beyond a value Qh of the 
charge, bubble dimensions may reduce to below the value 
of the van der Waals hard core radius for the gas enclosed, 
which is physically untenable. 

Hence, the value of Qh^ the physically feasible maxi- 
mal limit to the charge the bubble may carry, will be less 
than Qmax for a particular Rq. Moreover, it depends as 
well on the amplitude of the forcing pressure Pg, with 
Qh decreasing with increasing Pg and also with decreas- 
ing driving frequency. Figure ([7]) show plots of Qh as a 
function of Pg for three different driving frequencies, for 
Rq = 2fim. Below a certain value of P5, Qh becomes 



nearly independent of frequency as well as Pg . 

V. BIFURCATION DIAGRAMS 

That the driving frequency influences the bubble 
dynamics is unquestionable. Techniques of dynamical 
systems theory have been used for long in the litera- 
ture to understand bubble stabilit y u nder variation of 
parameters (see for example [19|, |2j427|. Parlitz, et 
al. [19| have, in their work, extensively investigated the 
frequency bifurcation diagrams for the bubble radius at 
various values of the driving amplitude pressure, Pg. 

In Figs.(|8lfTT]) we have shown the bifurcation diagrams 
for the maximum radial amplitude Rmax of the given 
time series of the bubble, with the driving frequency as 
the control parameter for Rq = 1.45/i77i, 2/i7n and bfim 
for uncharged and charged bubbles. 

The bifurcation diagrams for various sets of parameters 
are constructed by sampling the time series after making 
sure the transients have decayed, for every time period 
T = Xjv of the external acoustic driving pressure elapsed. 
These sample points are precisely the points of intersec- 
tion of the trajectories in phase plane with the Poincare 
cross section, and the orbit formed by the points repre- 
sents the Poincare map. The bifurcation diagram is then 
constructed by plotting the sampled points calculated for 
a range of values of the control parameter (frequency or 
charge) and then plotting it with the control parame- 
ter on the horizontal axis and the sampled points on the 
vertical axis. 

In the response curves where period-doubling bifur- 
cations occur, the branches always merge back to give 
period- 1 oscillations. 

The presence of non-zero charge on the bubble ad- 
vances period-doubling bifurcations with the driving fre- 
quency as the control parameter. This is demonstrated in 
the bifurcation diagrams in Figs. (8-9). For an uncharged 
bubble with ambient radius i?o of 2microns at a driving 
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FIG. 8: Bifurcation diagram of a bubble oscillator with driv- 
ing frequency as the control parameter, and for Ro = 2 /am 
and (top): Ps = 1.2Po (i.e., Ps < Ptr), and (bottom): 
P, = Pt^ = 1.3Po. Q = in (a) and (c), and Q = 0.2pC in 
(b) and (d). 



pressure of I.2P05 period doubling is first seen at around 
720kHz for the uncharged bubble, while the presence of 
0.2pC charge advances it to about 600 kHz (Fig.dH a,b)). 
We observe that there are no chaotic regimes present at 
least till driving frequencies of lOOOkHz for low driving 
pressures such as this. 

Figs.([8l fTQ|) show that increasing the external pressure 
Ps also has the effect of advancing the succession of 
period-doubling-period-halving bifurcations both for the 
charged as well as for the uncharged systems. For in- 
stance at I.3P0 (Figsl8]c,d), the first period doubling bi- 
furcation occurs at a forcing frequency of approximately 
320kHz, followed by period halving bifurcation at 350kHz 
leaving period 1 oscillations, whereas on introduction of 
charge Q = 0.2pC, the first period doubling bifurcation 
makes its appearance much earlier, at about 295 kHz, 
only to merge back to period 1 oscillations through a pe- 
riod halving bifurcation at 315 kHz. As one increases the 
driving frequency further, one observes the occurrence of 




200 400 600 800 

V (kHz) 




200 



400 
V (kHz) 



600 



800 



FIG. 9: Bifurcation diagram with respect to driving frequency 
of a bubble oscillator with Ro — 2/j.m and Ps — 1.4Po for 
an uncharged bubble (top) and having charge Q — 0.2pC 
(bottom). (Note that Ps > Ptr.) 

a sequence of period-doubling - period-halving bifurca- 
tions. 

It should be noted that the chaotic regions make their 
appearance at the upper transient threshold pressure Ptr 
(which for an uncharged bubble of Rq = 2/im is 1.3Po), 
and become more prominent for Pg > Ptr (Fig|9|). 
At large driving pressures, bands of chaotic regimes 
are present at high values of the forcing frequency, in 
agreement with observations of time series data. We 
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FIG. 10: Bifurcation diagram with respect to driving fre- 
quency of a bubble oscillator with Ro = 5/im and for Q — 
O.SpC. The driving pressure Ps = 1.4Po (Here, Ps > Ptr)- 
Successive images (b,c,d) show magnification of preceding in- 
sets showing period-doubling and chaos. 

show this in Figs.([9|) for Pg = 1.4Po • chaotic behaviour 
is seen to be present even at around 270-300kHz for 
charged and uncharged bubbles. 

In Figs. (p!Q|) the sequence of period-2 and period- 1 
oscillations generated is shown for a slightly larger 
bubble, with Rq = bfim and bearing charge Q = 0.8pC, 
driven at pressure amplitude Pg = IAPq. 

From these observations and other plots (not shown 
here) we deduce that the maximal radial amplitude 
of the bubble of a given equilibrium radius Rq shows 
chaotic behaviour as a function of the driving fre- 
quency v, for P. > Ptr at large values of i^. It 
was shown in J23| that Ptr = 1.1 2Po for Rq = bjj.m 
and Ptr = 1.3Po for Rq = 2fim for the uncharged bubble. 

For smaller frequencies, such as in the sonolumines- 
cent regime, the presence of charges do not appear to 
introduce period-doublings in the system. However, the 
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FIG. 11: Bifurcation diagram of a bubble oscillator with driv- 
ing frequency as the control parameter for a smaller bubble. 
Ro = 1.45/im and Pg = 1.4Po. 
(a) Q = 0; (b) Q = O.lpC. 



effect of charges in bringing about drastic changes in the 
bubble stability is more pronounced for smaller values of 
Rq. This is demonstrated in Figs. (11 a,b) for a bubble 
with Ro = 1.45/im and driving pressure amplitude Pg of 

1.4Po. 

The presence of O.lpC charge on the bubble (Figs. (fTTb )) 
induces a period-doubling bifurcation at a driving fre- 
quency of 880 kHz followed in quick succession by a 
period-halving bifurcation at 1040 kHz. These are ab- 
sent for an uncharged bubble (Fig. (fTTk )). 
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Tmax depending upon the amplitude of the driving pres- 
sure, Pg. 

At lower pressures, i.e, for Pg < Ptr (for example for 
Ps = l.lPo for Rq = 5/im), Tmax increases with driv- 
ing frequency. However, as this value of Ps falls in the 
vicinity of the transient threshold in the unstable regime 
(Fig. (4a)), we would expect Tmax to show large oscilla- 
tions with frequency, as is also seen in Fig. (13a). 

At higher pressures, i.e, for Ps > Ptr (for example for 
Ps = 1.25Po for Ro = 5/im), the maximal temperature's 
frequency dependence is the opposite, with T^ax decreas- 
ing uniformly with increasing frequency, showing oscil- 
latory behavior (Fig. (13b)). T^ax shows a frequency- 
dependence of the form 



FIG. 12: Bifurcation diagram of a bubble oscillator with 
charge Q as the control parameter, and with Rq — 2/jm and 
Ps — 1.4Po and at zy = 300kHz. Non-chaotic windows are 
seen, most prominent one centred around Q — 0.125pC 

In Fig. ([12]) we obtain the bifurcation diagram of a bub- 
ble of Ro = 2fim driven by sonic pressure amplitude of 
1.4Po and frequency 300 kHz with charge as the con- 
trol parameter. The choice of 300 kHz for the driving 
frequency has been made using the bifurcation diagram 
with frequency as the parameter (Fig. (9)) where the sys- 
tem is just beginning to get chaotic at this frequency. In 
Fig. (12) we see an interesting non-chaotic region centered 
around Q = 0.125pC with period-doubling and period- 
halving cascades. While constructing the bifurcation di- 
agrams care has been taken to ensure that we work only 
within the range of charges permissible for a bubble of a 
given ambient radius. 

VI. THE COLLAPSING BUBBLE: FREQUENCY 
&- CHARGE DEPENDENCE OF TEMPERATURE 

Investigating the maximum temperature as a function 
of the driving frequency, we obtain the interesting re- 
sult that there exist two distinct domains of behavior of 



Tmax = a4 X (i/ - tts) 



-4/5 



(22) 



a^ and a^ being constants with appropriate dimensions. 
This is understood by recalling that the temperature is 
obtained from 



3\2/3 



^<°'i'|le|) 



(23) 



for r = 5/3. Making the approximation that {H/Rq)^ < 
1 and also that {h/R)^ < 1 is sufficiently small at most 
values of P, we can approximate Eqn. (|23|) by 



T ^ T{0)Rl/R^. 



(24) 



Since at regimes at or near the Rayleigh collapse, R{t) 
cj^/^, it immediately follows that 



T{t) oc uj 



-4/5 



(25) 



Values of Tmax at lower pressures are less than that at 
higher pressures. Temperature Tmax rises steeply with 
pressure Ps after some critical value of the pressure that 
equals the upper transient pressure threshold. Temper- 
ature calculations have also been reported by Wu and 
Roberts [28| and Yasui [29| giving high magnitudes of 
the temperature. The presence of charge serves to reduce 
the pressure at which a maximum temperature Tmax is 
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FIG. 13: Maximum temperature Tmax vs. u for different 
values of charge Q at (a) Ps = l-lft (lower / medium 
pressure, Pg < Ptr) & (b) Ps — 1.25Po (high pressure, Ps > 
Ptr)'i the contrasting behavior above and below Ptr can be 
seen, (color online). 



reached. As seen in Fig. (5b), increasing the driving fre- 
quency shifts the curves to the right, i.e., the same maxi- 
mal bubble velocity Vmax that is obtained at some driving 
frequency Ui for a pressure amphtude Psi is reached for 
a higher frequency UJ2 > cji only at a higher pressure 
Ps2 > Psi' Figs. (5b) and (14) where Vmax vs. P^, and 
Ptr VS. Q and Tmax vs. Q plots respectively are shown 
for Rq = 5/im and 2/im for charged and uncharged bub- 
ble at different driving frequencies, illustrate the effect 
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FIG. 14: (a) Ptr decreases with increasing Q over all driving 
frequencies. At a given charge, higher driving frequencies 
imply higher Ptr- (b) The Ptr — Q curves when scaled by Ptr 
at Q = 0, Ptro, all fall on a master curve obeying eqn. ([26|) . 
(c) Tmax vs. Q curvcs plotted at pressure value equalling the 
corresponding upper transient threshold value Ptr for that 
particular value of Q, all fall nearly on the same curve. Plots 
are for Rq = 2/im. (color online). 



of driving frequency and charge on bubble velocity, tem- 
perature and the location of the upper transient pressure 
threshold. 

The net effect of charge is to raise the possible Tmax 
for a bubble in comparison to the uncharged bubble. A 
comparative plot of Tmax reached at the upper transient 
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threshold pressure Ptr as a function of charge Q is shown 
for a few values of the driving frequency u in Fig. (14). 
The temperature Tmax increases with charge Q for all 
forcing frequencies. 

The value of P^r (for a given bubble-charge, Q) increase 
with v. Ptr at a given driving frequency decreases with 
increasing Q. The dependence of Ptr on Q over all fre- 
quencies can be captured by a normalized plot of Ptr /Ptr o 
against Q, where Ptro is the upper transient threshold 
pressure value at zero charge, for a given frequency. This 
yields a master curve approximately obeying a relation 
of the form 



PtrO 



(26) 



'"^ 1 + 0.794Q2' 

as seen in Fig. (14b). The maximal value of temperature 
Tmax reached in a driven oscillating bubble, at the cor- 
responding, respective Ptr (which varies with u and Q), 
over all values of charge Q, seems almost independent of 
the driving frequency v^ as shown in Fig. (14c). 

At higher pressures, beyond the upper transient 
threshold pressure, bubble velocities become larger and 
of the order of c. We argue that the maximal kinetic 
energy ^MR^ ^ ^Mc^ would approximately equal the 
dominant electrostatic contribution to the potential en- 
ergy QV(47rei?o), so that using M = 47ri?gp/3, Q^ - R^. 
This argument is independent of the driving frequency at 
which the bubble is being forced. At a sufficiently high 
charge Qmaxi a bubble of ambient radius Rq will col- 
lapse to the same minimum radius Rmin-, independent of 
the driving frequency of the forcing pressure amplitude. 
However, this Qmax value will typically be greater than 
Qh^ the upper bound imposed on the charge Q by the 
physically realistic requirement that Rmin does not go 
below the van der Waals hard-core radius. Thus while 
this results in Qh being the greatest, physically realistic 
value of charge that a bubble can carry, we can still read 
off the value of a larger Q = Qmax from Rmin vs Q plots 
at high driving pressures, by identifying the Q at which 
frequency independence of the curves sets in and all the 



curves for different driving frequencies all converge to the 

same Rmin- More detailed discussions of Rmin dependen- 

i 

cies are included in our earlier work [23], and we do not 
show the plots here. 

Hence a comparative estimate of this maximal charge a 
bubble can carry, Qmaxi for two different values of the 
initial bubble radius Rq^ say, Rq^ and i?06, would be ob- 
tained from 



Rob 



(27) 



(for larger bubbles, of order 0(/im) and above). This is 
essentially a statement that the maximal surface charge 
density a bubble can carry is approximately same re- 
gardless of its initial ambient radius for micrometer and 
larger bubble radius, given that all other system pa- 
rameters like the surface tension, pressure conditions, 
viscosity, etc. remain unchanged, while the influence 
of the effective surface tension (including the correc- 
tion for charge present) is predominant in the submicron 
range. Indeed, for sub-micron and at very small bubble 
sizes where the surface tension and charge terms become 
just comparable, we would have instead (comparing elec- 
trostatic energy with surface tension or elastic energy) 
Q^/Ro ^ kRl = ttiujIRq^ with k = muol being an ef- 
fective spring constant and uo^ some natural frequency 
and m the oscillator mass, so that on substituting for 
m = (4/3)7rpi?Q, we get 

/ -^Oa \ ^ ( Qrnaxg \ /^on 

\R0cJ \QniaxcJ 

for two different values of the initial bubble radius Rq^ 
Roa and -Roc- 

That such a non-rigorous approach cannot give any ac- 
curate numbers is obvious. Nonetheless, it is useful in 
giving us rough estimates of the maximal charge that 
the bubble can carry, in the absence of a constraint such 
as that imposed by the van der Waals hard-core radius. 
A comparison of the numbers so obtained in this rough 
and ready way to that obtained from the numerical re- 
sults is given below in Table 1 for three different values 
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of ambient bubble radius Rq, at Pg = I.SSPq. 

TABLE I: Qmax ratio from Eqn.(27) (through Rq ratio) com- 
pared with numerical values of Qmax from data. 



Table 1 


Roa 


Rob 


ymacca 


Qmaxt 


i^y (fet) 


2 /im 


5 /im 


0.23 pC 


1.24 pC 


0.16 


0.18 


10 /xm 


2 /im 


4.7 pC 


0.23 pC 


25.0 


20.43 


10 /xm 


5 /im 


4.7 pC 


1.24 pC 


4.0 


3.79 



VII. CONCLUSIONS 

The presence of charge on a bubble suspended in a fluid 
influences the bubble's oscillations under ultrasonic forc- 
ing, and some of the aspects of the dynamics have been 
addressed in this work, taking the polytropic constant 
F = 5/3 which governs the equation of state for adiabatic 
heat transfer. A dimensionless constant ( which we intro- 
duced in an earlier work [23] helps us to identify clearly 
the Blake threshold and the upper transient threshold Ptr 
for acoustic cavitation. We use this to understand the 
influence of driving pressure Pg and frequency u of the 
applied ultrasonic field on the bubble oscillations. The 
presence of charge reduces the effective surface tension on 
the bubble walls so that its maximum radius Rmax at- 
tained during the expansion phase is larger than when it 
is uncharged; similarly the minimum bubble radius dur- 
ing collapse Rmin is much smaller in magnitude when the 
bubble is charged. The charged bubble undergoes a more 
violent collapse, achieving far higher temperatures in its 
interior in comparison with the uncharged one. We find 
that when Pg < Ptr^ the maximum temperatures Tmax 
achieved in the bubble increase with increasing v and 
charge. For Pg > Ptr-, Tmax obeys a power law decrease 
with respect to z^, with an exponent of -4/5. The power 
law behaviour is also obtained analytically through scal- 
ing arguments near the regime of Rayleigh collapse. 



Bifurcation diagrams of the maximal radial amplitude 
of the bubble as a function of the driving frequency show 
the presence of chaotic regimes for Pg > Ptr for any given 
ambient bubble radius at fairly large driving frequencies. 
The route to chaos is through period-doubling followed 
by period-halving bifurcations. The effect of charge is to 
always advance these bifurcations. At the lower end of 
the ultrasound spectral range, for instance in the sonolu- 
minescent regime, the presence of charges do not appear 
to induce any period-doublings. 

Consistent with the fact that the presence of charge has 
a greater dominating effect over surface tension on bub- 
bles of smaller equilibrium radii [23], the bifurcation di- 
agrams demonstrate that the effect of charges in dras- 
tically changing bubble stability is more pronounced for 
smaller bubbles. 

We obtain also the bifurcation diagram of the maxi- 
mal radial amplitude at any given Pg as a function of 
the charge at large driving frequency. Here too, period- 
doublings and period-halvings are seen interspersed with 
large chaotic regimes. 

We obtain analytically an estimate of the minimum 
charge Qmin required on a bubble at a given magnitude of 
applied pressure to attain a certain value ci of the bubble 
radial velocity. We find that this is related by a simple 
power law to the driving frequency of the acoustic wave. 
We show that above a critical frequency ujh^ uncharged 
bubbles necessarily have to oscillate at velocities below 
ci. The calculations are reproduced numerically also. 
Further, ujh depends upon Pg. 
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